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AN EXPERIMENTAL INVESTIGATION OF THE GAINS OF DISLECTRIC 
ANTENNAS 


INTRODUCTION 

The idea of using dielectric materials for the 
control of the direction of electromasnetic rediation 
was probably first conceived by Hertz in the Nine- 
teenth Century?, when he observed the strong influence 
of vitch blocks on this little understood form of 
emorey fT low. Subseattently, in the early part of the 
Twentieth Century a few investigations were revorted 
eoncerning the use of dielectrics as wave guides, 
but there were only occasional contributions to this 
theory. Rapid progress in this direction was not 
begun until the availability of microwave equipment 
in the 1930's. By 1939, dielectric wave guide theory 
had been thoroushly studied and a few experiments 
had been attemoted in the use of cielectrics as 
directional radiators, With the Beis of the war, 
and the ravidly increasing need for more highly 


directive antennas, dielectric rods were studied, 


lL. Superscrint numbers refer to references listen in 


Bibdlieszranhy. 
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tried, and put into use beginning in the early 1940's. 
In 19435 the U.S. Navy put into service a radar 
equinment, designed by the Bell Telephone Laboratories, 
which emoloyed an array of forty-two rods of 
Polystyrene to give a sharp directional beam in the 
horizontal plane. (Radar Equipment Mark VIII). 

The Germans, about the same time, began to use a 
directional antenna on aircraft’ which consisted of 
four dielectric rods arranged to sive a symmetrical 
pattern. Subsequently, a large number of dielectric 
antennas have been built, and an ever increasing 
demand for these small size, highly directional 
devices has encouraged more and more attention to 
their theory of operation, 

It should be pointed out in the beginning that 
the classical metnods of antenna design, based on 
Maxwell's equations, have not been successfully 
applied to dielectric antennas whose dimensions have 
been of the order of a wavelength. As clearly pointed 
out by Stratton’ in his book, "Electromagnetic Theory", 
Maxwell's equations are valid, in seneral, only to 
large scale phenomena, and attempts to extend them 
to the so called "Quasi-microscopic state" lead only 


to partial success. Although many of the laws of 





3. 
classical electrodynamics apply to submicroscopic 
domains, it is not possible to apply them generally. 
On the other side, there is. quantum electrodyna:rics, 
which at this writings, have not been extended 
through the transition areas. The statistical 
DCm aco GV Creo lota.c Nitmbher Ol vauzoms, of the 
equations of quantum electrodynamics ultimately 
must lead to the relations derived by Maxwell. 

Design of dielectric antennas to date has been 
based larzely on cut and try methods. Some of their 
characteristics have veen determined by various 
exverimenters, and considerable theoretical work 
has been done in an endeavor to explain their 
operation. But there is no mass of data for the 
theorist to draw uvon, nor is there an adecuate theory 
at the present time to enable calculation or prediction 
of the characteristics of a given dielectric shape. 

It is the nurpose of the writer, in undertakinz this 
work, to determine experimentelly the "concentrating 
action" of certain of these antennas, and to express 
028 in terms of sain relative to an isotronic 
oscillator; or, what amounts to the same thins, to 
express this gain as a ratio of effective area to 


actual area. It is believed that these data, thus 





4, 
obtained, may be of use in subsequent investigations, 
or for other investigators, to ultimately develop 
a satisfactory theory of the function of the dielectric 
antenna. Finally, available theoretical approaches 
will be discussed objectively in the light of the 


experimental results obtained herein. 


work of other investizators: 


Mallach®, in Germany, presented a thesis in 
1939 whicn revorted considerable work in the use of 
Gielectrics to produce a directional beam. He also 
produced a gimple theory eaten exovlains quabitatirebyy 
the action of the dislectric, but failed quantitatively 
tO verify all of the known data, Takines a dielectric 
filled wave guide, whose walls are considered 
perfectly conducting, excited at one end and open at 
the other, it is well Known that radiation is emitted 
from the open end, and this radiation is directional 
along the projected axis of the guide. If the 
conducting sheath is removed from a part of the guide, 
nearest the open end, leavinz some of the dielectric 
rod vrotruding, radiation will occur both from the 
end of the dielectric and from its side walls as 


well, Mallach's theory is based on the consideration 
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that each element of the surface acts as a radiator 
and the phase of tnis excitation devends on the 
extended length of the rod.. He assumes in this work 
that the wave is linearly volarized and of constant 
intensity and is transmitted down the rod without 
attenuation or reflection. Furthermore, he coes not 
consider the rod cross section area except for the 
determination of the wavelengtn, He also considers 
only small anzles of diffraction, takinz no account 
of the spherical nature of the wave. Mallach 
reported a considerable amount of exvnerimental data 
on both rods and dielectric tubes, and arrived at 
some simole formulas for choosing dimensions. For 
the sizes of rods and tubes thet he used, his simmrle 
theory and formulas are satisfactory, but are of no 
use in obtainins lobe pnatterns for the entire 
radiation field nor for computing gains. 

Southvorth!, at the Bell Telephone Laboratories, 
whose early publications on wave guides and wave 
guide theory are well known, has done considerable 
experimenting for more than a decade with tavered 
Gielectric rods, both rectangular and circular, of 
diameters of a helf wavelength or less. He and his 


associates have exverimentally established that 
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the gain decreases as the attenuation in the rod 
increases. 

WAlces?, in a thesis presented to the University 
of Paris in 1946, examined the problem from the 
point of view of the lens action of the dielectric 
on radiation of waveléngths comparable to the 
dimensions of the lens. He pointed out that the 
apparent index of refraction plays an important 
Bowe widetic enerey sain of the dlelectric raciator. 
By vlottingz experimentally obtained enersy gains 
versus lenzth of exnosed dielectric he has shown a 
correspvondence of periodic zain to’Snell's law of 
refraction. Furthermore he has shown that the 
apparent index of refraction is itself not constant, 
but denends on the exnosed length of dielectric and 
Pmeros> SoCtlon dimension. A formula for calculating 
the anparent index from the real index of refraction 
ig develoved in his paver. Wilkes also computed some 
radiation patterns based on analozy to more 
conventional antennas, which verified very nicely the 
experinental data for the major lobe. But, as his 
paver points out, his assumptions used would cause 


considerable errors in the side lobes. This 
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calculation assumes constant wavelength, which he 
shows to vary somewhat throushout the lens length. 
He also assumes, for this calculation, constant 
displacement currents. ‘the complexity of an exact 
formule tion of wavelensth and intensity would lead 
to expressions which would be extremely difficult 
if not impossible to solve. 


11,12, 16 have subsequently 


peveral investigators 
published attempts at calculation of radiation 
patterns. All of the recent attemots knovm to the 
author are based on the formal solution of the field 
equations given in Schelkunoff's well known book’, 
as aonlied to waves in dielectric wires, and exoressed 
in terms or cylinder functions. Elsasser in 1949 
published an analysis of the various non-radiative 
modes ina circular dielectric rod, in which he also 
pointed out a perturbation method of calculating 
dielectric loss. Blsasser, in this paper, did not 
attempt calculation of radiation from a finite length 
rod. This theoretical study is in substantial 
agreement with an experimental investigation of 
dielectric cylinders as wave guides performed by 


10 


Chandler at RCA laboratories. 
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Marston!+ of the U.S. Naval Research Laboratory 
has attemvted calculation of radiation from dielectric 
rods of uniform circular cross section. Ue employed 
the method of assumed fields based on Schelkunoff's 
Kauivalence eeoneren & published in the Bell System 
Technical Journal in 1936. Marston computed the 
radiation field arising from equivalent sources at 
the sides or the rod, and senarately for those at 
the ends of the rod, and points out that maximum 
Pow in the axial direction is achieved when the 
total radiation from the end is in phase with the 
radiation from sides (for modes of type n = 1 only). 
Marston's calculations also show that the actual 
values of gain for dielectric circular rods cannot 
in general be vredicted from ordinary array tneory, 
but ovrosvects of achieving this exist through the use 
of his complex eynressions for the components of 
electric field intensity. However, all of Marston's 
assumptions are not in asreement with findings of 
other investigators; for example he assumes constant 
wavelength within the dielectric body. If tnis 
variable were substituted for the constant, 1t would 


make his equations almost honelessly complex. 





A great deal of the knowm information on 
dielectric antennas has been assembled by We eneyae 
in @ dissertation rresented to the University of 
Texas in June, 1950. He renorted also a very 
substantial amount of experimental data, as well as 
@ slightly different method of attempted calculation 
of the fields. The method vresented is also based 
on Schelkunoff's Equivalence Theorems, but uses an 
extension devised by Horton, in which He comoutes 
fictitious electric and magnetic current sheets 
about the dielectric which are a useful aid in the 
analysis. Then he obtains the vector notentials at 
a point in space due to these sheets and commutes the 
E and H fields at the point. In his calculations he 
neglects the radiation from the end of the rod, 
neglects reflections in the rod, and neglects any 
Soanees Waiecm mieht occur to the wavelensth of the 
wave within the dielectric. McKinney's exnerimental 
data agreed with his theoretical calculations quite 
closely. For the sizes and shanes used and for the 
To mode, it annears that his solutions are auite 
precise. However, his renort does not show any 
calculations of patterns for dielectric antennas 


excited in the dominant mode, which is the mode most 
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commonly used, nor does his theory contain a precise 
statement of the optimum gains which may be obtained. 
The T™o1 mode sives a null on axis, which is 
unacceptable for most cases where dielectric antennas 
are reauired. 

| It should be pointed out that no theoretical 
treatment of any dielectric antennas other than those 
of circular cross section have been found in any of 
the literature. As will be brought out later in this 
essay, an attempt at a straightforward solution to 
the rectansular shaped antenna by any of the classical 
methods would lead to exceedinsly complex equations 
whose solutions would involve serious difficulties. 
Nevertheless, in a great many applications, the 
rectangular antenna is the one preferred due to the 
predominant use of rectangular wave guide which 
assures no rotation of the plane of polarization. 
Hence this essay includes experimental sains obtained 
for both rectansular and circular cross section 
guides, including such correlation to existing 


theories as are vossible. 





a, 
CHAPTER 1. 


REVIEW OF BASIC LAWS OF ELECTROMAGNETIC RADIATION 


The objective of Chapter 1 is to develop the 
ecuations for the provagation of electromasnetic 
maves in a circular dlelectric rod, infinitely 
lons, imbedded in a homogeneous medium. The starting 
noint will be Maxvell's equations. The development 
will be carried out for the cylindrical case because 
treatments of this are available in the literature. 
Throushout this work the MKS system of units will 
Ine used. 

It has been vointed out in the introduction 
that Maxwell's ecuations hold with certainty only 
for macroscopic vrenomena. Accordingly, we make 
the assumption throushovt Chapters 1 and 2 that only 
mMacrosconic conclitions are considered, that the 
media considered are isotrovic and homogeneous, 
and that, except for stated discontinuities, the 
ohysical properties in the neishborhood of every 
point are continuous. Under these conditions, the 
following Maxwell's equations have been firmly 


established. 


1) vxe +22:0 


12. 


Soe vty - 2-7 


where the vector quantities E, B, H, D and J are 
respectively the electric intensity, the magnetic 
induction, the masnetic intensity, the electric 
displacement and the current density. 

Since the divergence of the curl of any vector 
is identically zero, we obtain by taking the 


divergence of (1) 
eo) We 28 =O 

and under the condition of continuity imposed above 
i278 =0 


Since the divergence of the magnetic intensity does 
not change with time, and cortainly at some time 


it did not exist, it follows that 
(5) Wa Giee=.0 
Likewise, taking the divergence of (2) and 


commuting the partial and del signs 
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Now if we denote the charge density at any 
point by Pf, the scuation of continuity or the 
conservation of charge states that the time rate of 
change of Pnlus the divergence of the current density 


is identically zero, or 
eu 
Sie) te eae) 
Substituting for v-g in (6) 


(8) 2 V7D- =O 


or 
(9) 2.(V7-D-f)=0 


Applying again the argument that if the time rate 
of change of the quantity is zero, the quantity 
Ltself must be zero if it ever has vanished, it 


follows that 


Go y* D= PLP 


This tells us that the charges distributed with the 


density Pp constitute the sources of the vector D., 


meecuric and Marcnetic Polarization 
In order to describe the electromagnetic state 


of a substance, it is convenient to define two 


14, 


additional vectors, the electric and the magnetic 


polarization vectors. By defini:ion 


Om, Gf =)-€.F 


i 


Vf 


2) (e-em 


where € is the inductive capacity in free snace 

and 4, is the permeability of free space. Therefore 

in free snece, the volarization vectors vanisn. 
Substituting the values for D and H from (11) 

and (12) into (2) and (10), we obtain the following 


(13) VxB-4& $6 2a (y+ Sh +V7xXM) 





(14) V-E=¢(P-V-P) 


In any vhysical media which is isotropic and 
in which the field vectors are linear, the Hlectric 
Displacement and the Magnetic Induction are related 
to the corresvondings field intensity by Hes conaeene 


of the medium 


(15) D = €£ 
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(16) B=4H 
From the above equations (13-14) it is possible 
to account for material bodies in an electromesnetic 
field by an equivalent charge of density -VP and 
an equivalent distribution of current density 23+ 7*M 
It is convenient at this point to define two 
edditional coefficients of pronortionality. The 
electric and magnetic susceptibilities of any media, 
excluding ferromagnetic material, will be denoted 


BX 5 andx,, and are defined by 
CMe = Xe bak 
Oe ay, 


We can define these scalar quantities in this way, 
since in isotropic media, the polarization vectors 
are parallel and provortional to the corresponding 
field vectors. 

It 1s also convenient to define the rolative 


exselectric constant 
_ € 
(19) f = eo 
and the relative permeability 


(20) Kn = 








Sy 


From (11), (15), (17), and (19) 4t follows that 
fen 2, = € =i 
and from (12), (16), (18), and (20) that 


Ca, Xt, =i! 


Vector and Scalar Potentials 

To facilitate the analysis of electromagnetic 
fields it is useful to define some auxiliary functions 
called potentials. These may be shown as follows: 

Equation (5) shows that the divergence of the 
Magnetic Induction Vector, B, is zero at all times. 
According to this, the field of Bis always solenoidal, 
and therefore B can be represented as the curl of 


another vector A. 
(23) B= VXA 
Substituting this value of B into equation (1) 
gives 
2 
(24) Vx (E+ $4) = 0 


Thus the field of the vector (E+ 24 ) is 


irrotational and equal to the gradient of a scalar 





ae ts 
function which we may call p. 


> 3A 
(25) +y24=-v¢ or B=-V9-3¢ 


The function A is called the vector potential of 
the field, and the $ is the scalar potential, 

It should be noted, however, that the vector A 
and the scalar p are not unique. Thus B could have 


been represented by the curl of another vector 
(26) B= VXA 


where 


i} 


(27) A=A,-VY 


where Y 1s any arbitrary scalar function of position, 


Lt would similarly follow that 
(28) E=-V¢, - 2% 

and @ and }, would be related by 
(0) =, + ae 


Likewise an infinite number of potentials describing 

the same field could be derived in ths same manner. 
By the definition of &, equationg (15), and 

using the ahove defined potential, equation (25), 





eS. 


Es A 
(30) ee en 


And by ecuations (16) and (23) 
(31) H= ~ZV*A 


Substituting (30) and (31) into eouation (2) 


vx(kvia)- 3f-e(vd +34) ] = 


and multiplying the above by “ and commuting 2, 


with V gives 


od : 
(32) Vai +uev se tue rh ee] 
By substituting (30) into equation (10) and dividing ° 
py -é 


(33) vivg+ v 242-2 P 


Now all varticular solutions of (32) and (33) 
te2d to the same electromagnetic field when subjected 
to the same boundary conditions. They differ among 
themselves only by the arbitrary function Y . In order 
to synecify the functions more precisely, the 


condition may be imposed that A and g also satisfy 





We, 


(34) WA +e 2o =O 


In order to do this, it is necessary that W satisfy 
2 = 
(35) V*W - ue Eh = VA. tuk 2h 


Where A. and %, are varticular solutions of (32) 


and (33). 
The vnotentials A and p by the aid of the above 


are completely svecified., They are solutions of 
the following ecuations, which are the same as (32) 
and (33) excent for the substitution of VA =-ue 2h 
from (34). 


(36) VX VKA - VV-A +48 24, = AT 
(37) VW°¢ -ubé 28 =-7P 
By use of the vector identity 


(38) VWxXVXA = VVA-VVA 


ecuation (36) reduces to 


(39) VWA-xe 24 =-ay 





20. 
The relations (23) and (25) for the vectors 


Bend = are not senerel. To them may be added any 


morulLecular SOlUtTLONM Of the homogeneous ecuations? 


(40) WxE + $2 = 0 


(41) VXH - 3 = O 


(42) V' B= 0 


(43) VY: D=0 
From the symmetry of this system +t is evident that 
it can be satisfied identically by 


(44) D= -VxA" 


¥ 
(45) H= -vér-S6 
) 
from which may be obtained 


(46) m= -~ 49x 


( 47) Bs -U(Vg*+ 24) 








' Pls 


But the new cotentials,A* and *, must be subject to 


the conditions 


ae 
(48) viAr-x“e af = 0 


& 
(9) gt ae Bf 


i 


(50) TV At tue 2& 


Tt 
© 


A general solution of the inhomogeneous system 


given by equations (1), (2), (5) and (10) is, then 


(51) B= VYxA ~4 2-4 vgr 
(52) Es-vg - 4- evra" 


The functions A” and g* are potentials of a source 
distribution entirely external to the region 
considered. Usually An and >” are set equal to 

zero, and the potentials of all charges, wherever the 
source, are represented by A and $. 


At any voint where the charge and current 
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densities are zero, a vossible combination of 
is 
potentialsAA, = O and 5 = Og», The function is! 
then according to (35) any solution of the homozeneous 


eauation 
(53) V*y - ub =H =0 


Since $, = 0 satisfies this same ecuation at this 
point,Y¥ may be chosen in such a way that p vanishes. 
In this case the field can be expressed in terms of 


only the vector potential. Equations (23) and (25) 


become 
(54) Be VIKA 
aA wv 
(55) B= 7 oF 


and (39) and (34) become 


a 


(56) VA-ué 4 =0 


(57) V'A =O 


The foregoing varagraphs have shown how an electro- 
magnetic field in a homogeneous medium may be 
reduced to a vector and a scalar potential which 


will satisfy Maxwell's equations. 
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Hertz Vectors or Polarization Potentials 
Lt is vossidle under ordinary conditions to 
define an electromagnetic field in a homoseneous 
medium by means of a vector function alone. 
Considering at first, no free charges or 
conduction currents within the medium, the field 
ecuations reduce to the homogeneous form (equations 
(40) through (43) ). Assume 2 vector, II, whose 


time derivative is yvropvortional to the vector 


potential A 


(58) A= ué = 


and substitute in (23) and (25) 


(59) Be Ue Vx et 
(60) es -VO-ué Sh 


Using the relations B =uffand D =&£,substitute the 
above values in the homogeneous equation (41). After 


commuting the sy andy signs, the result is 


(61) 2 (Vx Wx + V6 +u€24) = O 


4 
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Since the scalar function d is arbitrary as long 
as it satisfies equation (34), it may now be chosen 


to also satisfy 
(62) @=-VI 


if the above cnosen value of d is substituted in 
equation (51) and the result integrated with respect 


to time, the following is obtained 


aa 
(63) VxVxZT -vvz tuk SF: =a constant 
¥ 


Toe value of the constant does not affect the 
determination of the field, and so it may be replaced 


by zero. very solution of the vector equation 
az =f 
(64) VxaVal - VV tué sm 


determineé an electric field through 


(65) Bw Vi 
and 
(66) Es VVI-ué SH 


Equation (65) is identical to (59) and (66) is the 
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same as (60) with the substitution of YY 7 for 
its equal -Vo% by virtue of (62). 
Applying the vector identity, ecuation (38), to 


ecouation (64) gives 


(67) Vi VE -Me =a -, = 0 


The requirement that p satisfy equation (53) is 
met in ecouation (64) and (67). 

Inasmuch as the Electric Disvlacement Vector, 
D, as well as the Magnetic Induction is solenoidal 
in a charge free rezion, an alternative solution 
could be similarly constructed of the form 


en * 
(68) at = wé 2% 


ik 


oO 
ts 
i} 


(69) 


¥ 
(70) D=- 46 Vx st 


VV 0" ~ 4b 235 


an 
it 


(71) 


where II is any solution of (64) or (67). 
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It therefore follows that the electromagnetic 
field throughout a homogeneous isotropic region in 
woich no charge ior current exist could be resolved 
into two partial fields, one derived from the 
vector II and the other from II”. The source of 
these fields is exterior to the region. In order to 
show the vhysical significance of the Hertz vectors 
it is necessary to relate them to thelr sources. 
This may be done by finding inhomogeneous ecuations 
from which (67) is derived. 

The vector D may be expressed in terms of the 
electric field intensity and the electric polarization 
vectors, using equation (11), as D =@6, E +P. Then 
substitute this in the homogeneous field equation. 


Eouation (41) becomes 


(72) V«H -& 2% = 2 


; % 
and (43) becomes 


(73) WE =-#& TP 


The above two equations as well es (40) and (42) 
are still identically satisfied by (65) and (66), 
provided@é is renlaced by a and that II is now any 
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solution of 
or 
Wee Vo Vi -~4b'o sae =e FP 


The source of the vector II and the electromagnetic 
field which it prescribes is a distribution of 
electric polarization P. 
The vector P is the electric dipole moment per unit 
volume of the medium. The vector II is associated 
with 2 distribution of electric dipoles and may be 
called the electric polarization potential. 
Similarly the field associated with II is 
established by a distribution of magnetic polarization. 
The magnetic induction vector, B, 1s related to the 
Meenetic field intensity, H, according to equationé 
(12), by B= 4% (H+ ™M). Substituting this into 
(40) and (42) gives 


wae \/*X Ee + u, 24 = 


x 
k 


and 


oa Ven = - V/V 


If it is further stipulated that II” shall be a 


solution of 
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GMs yay i" -u,6 28 =-/ 


and that 4. replace Xin equations (70) and (71), 
then eauations (75), (76), (41) and (43) are 
identically satisfied by (70) and (71). 

The electric polarization P may be induced in 
the dielectric by the field E, but it may also 
contain a part whose magnitude is governed by 
entirely external conditions. Insofar as dielectric 
antennas are concerned, only the indenendent vart, Po» 
4s of interest, which represents the electric moment 
of divole oscillation activated by external power 
sources. This 1s also true of magnetic polarization. 
To represent these conditions, equations (11) and (12), 


expressed in modified form are 


(78) E= ¢(D-f,) 


zB - Me 


v 


(79) H 


in which - and MS are vréscribed and indenendent 
of E and H and where the induced polarizations of the 
medium have been absorbed into the parameters & and 
A. The electromagnetic field due to these distri- 


butions of P, and M, are given by 





ao). 


27 * 
(80) B= VV I-uwé 2s -avn 8 
(81) Ha EV KS Z+VVI*-a £2 


Equation (80) above was formed, by substituting in 
ecuation (78) the value of D ziven by (70) and ES 
given by (74), replacing 4, and &, by “ and @. 
Equation (81) was formed by substituting into (79) 
the value of B given by (65) and M, from (77). 
For these ecuations to be valid, however, it must 
also be stivulated that II and II” be solutions of 


equations corresponding to (74) 


(82) VeVi ee oF =-4P, 





(83) VV: Vi." ie sn =~, 


Using the vector identity, equation (38), and 
substituting (83), it 4&%salso possible to write (81) as 





BOs 


(84) H = EV x2F + Vx Vx" -M, 


The last relation shows that, since B = VXA, the 
vector potential A can be derived from the Hertzian 


vectors by setting 
(85) Ae wh ruVrd—*-Vy 


The associated scalar notential function is 


(96) ge - V+ 


in which must also satisfy 





(87) Vp -ué S4 =0 


Boundary Conditions - Discontinuities in the Field Vectors 
The above ecuations have all been based upon 

the assumption of continuous and isotropic media. 

It is now necessary to consicer the relations which 

arise due to the boundaries between two media, where 

there occur abrupt differences in the inductive 

canecities,uandé, and in the conductivity, 6 . 

| To aid in the analysis, it is useful to postulate 


avery thin transition layer, 5, between the media, 
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wherein the parameters“4,é, ands, vary rapidly but 
continuously from medium 1 to medium 2. Within this 
layer, 5, the field vectors and their first derivatives 
are continuous, bounded functions of time and svace. 
Consider a small right cylinder, Figure 1, in the 
layer such that its end surfaces lie in the surfaces 
of the layer and sepnarated by a length Ag, 

From the volume integral of ecsuvation (5) and 
the divergence theorem of vector analysis’? it 


follows that 


(88) [VB dir = J B-n da =a 


where the second integral is a surface integral 

over the walls and the ends of the cylinder, hn being 
the unit vector normal to $ . If the elemental 

area of the end of the cylinder is sufficiently small, 
B may be considered to have a constant value over 

the end. Neglecting differentials of higher order, 


(88) may be apvroximated by 


(89) (B-n, + B-n,) A@ + contribution of walls = 0 
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(2) fed, G2 
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where 4. is unit vector normal to the surface in 
the sense inward to medium 1, and h, the unit 
vector normal in the sense inward to medium 2. 
bemec the contribution of the walls is vroportional 
to Ak , it will become vanishingly small when the 
limit is taken Ak 0 , and the transition layer 
shrinks to the surface 5S. Designating the value 

of magnetic induction at the point on 5 in medium 1 
as Bas and in medium 2 as Bas 
limit AL>O and AQ-O 


and vassing to the 


fo (BB) = O 


Thus the transition of the normal component of B 
across any surface of discontinuity in the medium 
feo) CONtinuows, 

Applying the same process of reasonins to the 
electric displacement vector D, equation (10) is 
integrated over a volume and its ecauivalent from 
the divergence theorem substituted, In this case, 
however, since V:D 1s eaual to the charge density, P , 
when it is integrated over'a volume it becomes 


equal to the total charge within the volume, aq, thus 


(0) [Dadar = [D-nda = 








Bor 


This charge,a, is distributed throughout the 
transition layer with a density p, and as the 
cylinder is reduced toward the limit, the q must 
remain the same because of the conservation of 


enarge. Thus 
(92) = Pal Aa 


As ~aAk7O, f must then, It is then convenient 
to replace the product, Fat » by a surface density 
Or charse per unit area,w. The trensition of the 
normal component of the electric displacement vector, 


D, across any surface, 5, is then given by 


(93) (D,-D)n = 


The presence of a layer of charge on 5 results in 
an abrupt change in the normal comnonent of D, the 
amount of the discontinuity being equal to the 
surface density in coulombs per squere meter. 

It is also necessary to investigate the 
behavior of the tangential components of E and H at 
the boundary. To aid this anelysis, consider a 
rectangle of length AS and thickness AX, instead 
of a cylinder, lying at the surfaces of the two 
adjoining media, Let 5, be the area of this 
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rectangle and C, be its closed contour. 
Takinz eouation (1) and intesrating over the 


surface 


(oa) J (FE) nda +[ nda =0 


where f, 1s the unit normal vector in the conventional 
positive sense, es shown in Fisure 2. By Stokes 
theorem”, this first surface integral is equal to 

a iline integral of the vector E about a closed 
contour enclosing the surface, Thus (94) may be 


expressed as 


(95) [Ends +f SF on, site f) 


Designate by 7 and fT, the unit vectors in the 
direction of circulation as shown in Figure 2. Then, 
neglecting higher order differentials, equation (95) 


may be expressed by 


(96) (£ | +E%) 4S + contribution from ends=- 2% 7, asa 








NG ane 
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But the contributions from the ends are proportional 
to AX which becomes vanishingly small when the limit 
is taken. Let n be the unit vector normal to 5 


as shown, and 7T be the unit tangent vector such that 
(97) (ena 

Also, it is a vector identity that 
(98) woxn-E =n nee 

Thus, when the limit is taken Af>o , AS-—o0 


(99) -[ne(E,-E,) + Lum (38 48)} =0 


Since f,#0 and the orientation of the rectanzle 
1s arbitrary, the bracketed exoression above must 


be zero, and so 


OB 
(100) Nn 4 ~ }=-tim Bo 
(E, Sy Lim ot X 
And since we have imvnosed the conditions that the 
field vectors and their derivatives be bounded, the 
right side of the above equation vanishes when AJ~9Q 


Therefore 


(101) nx«(E,-E£,) =O 





io 


and it is seen tnat the tansential comoonents of 
the vector # through a surface of discontinuity 
are continuous. 

The behavior of the magnetic intensity vector, 
H, at the coundary can be deduced by the same 
process. Intesrating eauetion (2) over the surface 


of the rectangle described above and avplyinzg Stokes 


theorem to the first term, there is obtained 
pie) a= Se? aga = ff S-naa 
eS aaa i ° 
and it similarly follows that 
: aD 
(103) nx (H, -H,) = Lum ($3 +J) Ag 


The limit of gf AX as A> oO will vanish 
because D and its derivatives ere bounded. But it is 


possible that Lum JAX may not vanish, because 
6 


the current I = J-nasakis squeezed into an 
infinitesimal layer on the surface as tne sides 

are brousht tosether, and J~> © , It is convenient 
to revresent this surface current density by K, 
which 1s defined as the limit of Jak as AX 


aporoaches zero and J becomes infinite, Then 





Dies 


(104) na(H,-H) =K 


When the conductivities of the adjoining media 
are finite, no surface current can exist, because 
E is bounded and the vroduct of 6fafl vanishes 
with af, In the analysis of dielectric antennas, 
this condition exists, and so (104) can be reduced 


ZO 
(105) nx(H,-H,) =0 


In the foregoins, expressions were developed 
for the transition of the normal components of B 
and D, equations (90) and (93), and for the 
tangential components of E and (H), eauations (101) 
and (105). By using (15) and (16), the normal 
components of H and E are immediately obtained for 
the transition of an electromagnetic field from 


one medium to another. 


(106) #(H,- 4H) =O 


(107) n-(E, -#E)= & 
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and for the tangential components of D and B 


(108) nx(D, - $ D) =0 


n«(B, ~ 9 8) =0 
(109) «(I 2 A, ) (for finite 6) 


Poynting's Theorem in rezions free of Hysteresis 


Poynting's theorem-gives the relation between 
the rate of change of energy stored in a field and 
the energy flow. This may be obtained as a general 
integral of the field equations (1), (2), (5), and 
(10). Sinee —-J has the dimensions of vower 
expended per ee volume, 1t is anparent that this 
may be obtained by a scalar multiplication of (2) by 


RB. Thus 
(110) E:VxH - £32 =E-S 


Likewise eauation (1) can be siven the units of 


power ver unit volume by multivlying it by H, Thus 


(111) H-VxE +H 3 =0 


ee ee — -—= eee 





bo. 
Subtracting the above two equations 
3D OB _ 
cue Se Vat se tee —- bw 


But, by vector identity, the first two terms of 
(112) are equal to Y-(ExH)and so this equation may 


be written as 
(113) VW(EXH) ted = -ES2- 12 


Integrating over a volume V bounded by a surface 
Db, and applying the divergence theorem to the first 


term, this becomes 
(aay [ (Exn)-nda + fe dav = -[(E-32 +H 8)ae 


The terms on the right above represent the rate of 
decrease of the electric and magnetic energy stored 
in the volume, The second term on the left expresses 
the sum of the losses in Joule heat and the power 
exnended by the flow of charge against the impressed 
forces. If there are no transformations of 
electromagnetic energy into elastic enerzy of a 


stressed medium, the balance is maintained by a flow 








of electromagnetic energy across the surface bounding 
the volume. This is Poynting's interpretation of 

the significance of the surface integral mr 4 214) 
The reduction of electromagnetic energy stored in 

the volume is partly accounted for by the Joule 

heat loss, partly by the energy introduced by the 
impressed forces, and the remainder flows outward 


across the bouncing surface 5, which is represented by 


(1s) [ S:nda = f (ExH):nda 


Now the Poynting Vector is defined by 
(116) S=EXxH vatts/meter* 


which may be intervreted as the intensity of energy 
flow at a point in the field, or in other words, 
the energy ver unit time crossing a unit area whose 
normal is orlented in the direction of the vector 
ExH . It is Poynting's theorem that the vector 
product (116) at any point is a measure of the rate 
of energy flow per unit area at that point. 

Another useful form of Poynting's theorem is 
the differential form, in which the Joule heat 


dissipation is separated from the nover exnended 
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by the flow of charge azainst the imoressed forces. 
While E reoresents the total electric ficld 

¢ 
intensity, let ©& be the intensity of the impressed 


electromotive forces, Then 


i 


(117) Fe -E 


J * | 
(118) See = 
Substituting (118) into (113) and replacing D by 
its equaléfand B bywH, the differential form 


becomes 


(119) 
W(ExH) +s? + S(Eet+ Aue) = Clg 


The Complex Poyntinz Vector 

Let h =MtW revresent the density of the 
electromagnetic energy at any instant and QotI°-E J 
be the nower expended ver unit volume in Joule heat 
and by the flow of charge against impressed forces, 
From (119) with these substitutions, the Poynting 


theorem from (113) can be written 
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(1200) oS + +O =0 


In a stationary field, h is ner endent of time and 
the middle term would be absent. Moreover, when 
the field is a simole harmonic variation of time, 
the time average of the enerzy density over a pveriod 
is constant and 2h = 9 , the bar over the symbol 
indicating an averase value for a period. In the 


case of periodic fields, it follows that 


(121) or 
J5-nda +f { Par ae) 


Now these field vectors are harmonic functions of 
time and can be represented in comnolex form with 
real and imaginary parts. By complex variable 
theory, it can be shown that the time average of 
the product of the real parts of two commlex 
quantities is equal to one half the real part of 
the product of the first quantity times the 


conjugate of the second quantity. symbolically 


(122) lA) AA.) =+ B(AA,) 








{ac 


where A, is the comolex conjugate of Ay. 
According to the above, the mean intensity of 
tne energy flow in a harmonic electromagnetic field 


is the real vart of the complex vector 4 Exe 


S = GlE)KG UH) 
(123) ie 
5 = + G(ExH) 


This is the complex Poynting vector which may also 
be denoted by 


(124) sé stEx} 


In a plane wave the magnitude of H is + ,» where 
nh is the impedance of the medium, Therefore in 


engineering vractice the expression 


ef 


is frequently used, corresponding to the vector 


Coie 


(125) Power + 


expression (124), 








Representation of a cylindrical field by Hertz Vectors 











In a homogeneous medium which is free from 
any fixed polarization, the Maxwell Equations (1) 
and (2) are satisfied by the following relations in 


terms of Hertz vectors 


leo) mee yy X VL ay VX Ls 


fi 


(127) H = Valé +e) + Vx Vx" 


Mie first of these can be readily verified by taking 
equation (80), applying the vector identity (38) 


to the first term, and noting that ails —U gz =9 


by (67). The second one is similarly verified. 
Now consider a set of coordinate surfaces, 

Figure 3, formed by a family of cylinders whose 

elements are parallel to the # axis. These 

cylindrical surfaces need not necessarily be 

1? i. » and 1, 

are assigned as shown, 1, being normal to the 


circular or closed, Unit vectors i 


Zz 
cylinder, is tangent to it and parallel to theraxis, 


and Lo 18 Valens to the surface and perpendicular 








CIGURE 3, 





Relation of Unit Vectors to a Cylindrical Surface. 
The scenerators of the surface are varallel to the 


7 axis. 
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bo) 2 and 1. , forming @ risht handed orthogonal 
1 3 

set. Position with respect to axes determined by 

the three vectors is measured by the coordinates 


Be, , £ . An infinitesimal line element is 


(128) ds=i,h, dul +if,du? +i, dz 


Consider an electromagnetic field in an isotropic, 
homogeneous medium, undounded in extent, existing 
such that the Hertz vector, iI, is directed along 


the Z axis. II, =O II, = 0, and Il, =W, 


II = 0, The electric field intensity vector of 
this field, as ziven by (126) is 


(veo) EY = YXVKI 


and the magnetic vector of the field, by (127) is 
(130) H° =(€& +6) Vx 


The comoonents of the vectors E can be calculated 
from (129) by means of the vector formula for the 
curl in an orthogonal system of curvilinear 
coordinates (reference 19, pase 90 ), apvolying this 


twice. 








(131) 
HK Fi) - tele) pita (tele — 


OKs Fy fis Fa) 
“ETL \4 fs (hf 74 3 DL | ye fa aK, 


Bo 
Ar Ark; oa 3 
Where A, = gh Ao = a = Se 


this comes out to be, efter 


eeace id, -[,=O , 


taking the curl twice 





(izes ul Tz 
EY et JE at | 
o . | oe 


(152) z ~ A OF aye 


ef = plan B38) + G9) 


to compute the components of H from (130) 


Likewise, 
Simply take the curl once and multiply by é(% +) 


thus 











ofa 


) pot 


0 
He = (ez +6) 9 


(Ga Af = - (e2 +e), Ue 


Hf) = 0 


From tne z scalar comnonent of a function 
I,-’, an electromagnetic field has been developed 
which is characterized by the absence of any axial 
or longitudinal component of the magnetic vector. 
It was earlier nointed out thet II is the electric 
polarization potential, and therefore this field 
may apnrovriately be called a field of the electric 
type, but following the terminolosy of Schelkunoff’, 
it is more often referred to as the "transverse 
magnetic field.” 

It is necessary to digress to develop here 


the wave equation, which the function J[[, =¥ 


Must satisfy. 
Taking the curl of equation (1) and substituting 


AH soy als! 


(134) 9 Wx(VXE) = -4 (XA) 
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Since J =6f and D=€E , equation (2) can be 


written 


’ Substitute (135) into (134) and making use of the 
vector identity (40) 


oF 
(136) v(vE) - VVE= We SR ue dé 


Rewriting equation (10), substituting ¢£ for D 
(137) VE =e 


since in free space or in any conducting medium ? 
is independent of the field distribution, the first 
term of (136) may be taken to be ecual to zero. 


Then (136) becomes 


(138) TE -HE xe -AG a =O 


Since -J=6£ , this could also be written 


(139) VI ese, 46 2 =0 
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By eliminating E from Maxwell's equations, there is 


also obtained 


24 , . 
(140) V*4 -aé Sax —ug 24 =o 
Now, returning to the eauations of the cylindrical 


field, 7,=¥/, the components of E and H given in 
(132) and (133) must satisfy the scaler wave enuation 


(141) vy -4é 2, 46 SE =O. 


The Laplacian of the scalarwy, expressed in general 
orthosonal coordinates (reference 19, page /50 ),1is 


(142) 


VY= FF al a &, = aut ae Ge ace) + sie He oe) 


Substituting this in (141), where Az ee 
(143) 


a3 =, (& #) + zo 5 hs  %,) ise ~4ent ee 0 


The elementary harmonic solutions of (143) 


are of the form 
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ma t hz -AWX 
(144) Vy = Lh ee) e AKHE 


where f (uu?) is a solution of 
(145) 


| 0 (#2 of | 0 (& af 
ky hz (? 34) T PR au? & 2) + [R’-A}f = 0 


Similarly, a partial field can be derived from a 


a 
second Hertz vector II by the onerations 


(46) FU =-4% vx" 
A? = vxoxn* 


If the vector II is directed along the ¥ axis, 


the vector components come out to be 





(Caer 
E, - Aan a~pad? 
(2) WH Sie 
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ae 
(148) H,” = ore oe 








() | | 2 emda Os | 4: al) 
fi, = Af, - 9u"\ FR, ons) + aes ots | 


The scalar function irs Toeseoomocion of (141) 
and is a field of the magnetic type or transverse 
electric, characterized by a zero lonzitudinal 
component of k, 

If the partial fields derived from II, and IT, 
are suvernosed, the electromagnetic fields are of 
such generality that they can satisfy any vrescribed 
set of boundary conditions on any cylindrical surface 
moose senerating elements are parallel to the Zz 
axis, In general, the choice of the coordinates is 
limited pnrectically to systems in which the variables 
of equation (145) are separable. It is for this 
reason that thts development of equations of electro- 
magnetic waves in a rod is limited to the case of 


the circular cylinder. 





Dee 


The simvlest of the seoaradle cases is that in 
which the family y=constant is represented by a set 
of coaxial cylinders, where usr, u’=6 , &A,=/ 
and #,=F . See Figure 4, In this case equation 
(145) reduces to: 

(149) , 
(r+ + @-AF=0 
which can be separated by writing f(r @)as a vroduct 





(150) f =£(r) file) 


where f(y)is an arbitrary solution of the ecuation 
(151) 
tae (rE) H(A Aral = 0 

and f, (6) 1s an arbitrary solution of 


(152) qh aaa o 





Solutions of (151) are Bessel's functions, and 


LMG 


@ solution to (152) is (98) =e . The parameter, n, 
is a separation constant. Its choice is governed by 
the nhysical requirement that the field must be 

single valued at any fixed point in space. Since at 


the vresent, the discussion is limited to continuous 
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ates of the Circular Cylinder 
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media, the field is pneriodic in © and the values of 
Poerccmiievwed CO 1uverers. Of 1). 2, 5..5..) ClLther 
mOsSitive or negative. However, ina field represented 
by perticular solutionsof (151) and (152), which 
is bounded by planes 9=@, and @=@G@ , non-integral 
values would in general have to be assigned to n. 

abe weonstant h ls, im general, a complex quantity 
mervainin= to the propagation of the wave and is 
called She propagation factor. An explicit expression 
ror h in terms of the freaquency and the constants 
of the medium can be obtained only after ‘he behavior 
of ¥ over a cylinder f=constant on a plane Z= 
constant has been prescribed. 

The equation (151) is Bessel's equation, the 
Solutions of which are the Bessel functions v7, (VRUR r). 
In eccordance with terminology emoloyed by Jahnke 
and Emde and several other authors the term "circular 
cylinder function" will be used to denote any 
particular solution of (151) and desisnated by the 
letters f= Z, (@# 1) . The order of the functions 
whose argument is (Pg? r isn. Particular solutions 
of the scalar wave equation (141) which are periodic 
in t and @ can be constructed from elementary waves 


of the form 
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(153) 


De e SMEG oe (RB r} atthe ~4{Wt 


Reference is made to standard texts in 
mathematics for the pnroverties of the various Bessel 
functions, but it is essential to this develonment 
to note several pertinent characteristics. The 
Bessel function, or circular cylinder function of 
the first kind, is finite and converges for all 
finite values of the argument. For non-integral 
values of n it has two linearly inderendent 
solutions for J, (f) and d,(P) 3; but if n is an 
integer the two solutions are not independent. 

A second indevendent solution can be obtained by 
the Bessel function of the second xind, N., ($) , 
defined by 


(154) N,.(°) = sean LJ, Plcosnir ~-L,,(P)] 


When nis an integer, N_(f) as a second solution, 
but this is useless at On for the function becomes 


infinite at this point. Solutions can be obtained 
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in terms of still other functions, the Hankel 

(i | 
hme bons, fle (P) and H” (¢) , which are merely 
linear combinations of the first and second Bessel 


functions as follows 


(Se) Hy (P) = d, &) +iN, (0) 


(156) Hf) = JP) -iN, (P) 


But the Hankel functions likewise become infinite 
for f=0 and so are excluded from use where the 
argument is very small. #xpansion of the J and N 
functions shows a similarity to the cosine and sine 
functions and, as such, represent standinz waves 
when used with the time factor e“’*% . The H 
functions correspond to exponential functions with 
imasinary exvonents, and, when used with the time 
wt 


factor e represent trevoeling cylindrical Weves. 


The Field of Gircularly Cylindrical Wave Functions 


‘See Sa 





Every electromagnetic field within a homogeneous 


isotropic domain can be revresented by ‘linear 


Do: 


combinations of the following functions 


(157) Vn RA os pee (VRE re +t iRZ~Aiwk 


(158) VA = ail H (EP i) ptthz-iwk 


The first equation, (157), alone, is sufficient to 
represent an electromagnetic field within finite 
domains including the axis f=9 . But at great © 
distances from the source, (158) must be used, for 
as (WPK fl) becomes very large, the function reduces 
to a wave travelling radially outward. This can be 
more readily understood when it is noted that for 
very large arguments (f) the Hankel function is 


aporoximately given by 


(159) 
| i(p ~ 2th 
10) = Vee (ORD 


Bach elementary wave is identified by the parameters 


n, h, kK. When n is zero the field is independent of ® 





of. 


and is symmetric about the axis, when h is zero the 
propasation is vurely radial and the field becomes 
only two dimensional. The two functions, (157) and 
(158) may represent inhomogeneous plane waves as 
well as homogeneous. The nlanes of constant nhase 
are propazated alongs the # axis with a velocity 
Y= where « is the real part of h, but the 
amplitudes over these nlanes are functions of & 

and 6. 

Tne formulas for the comronents of the field 
vectors can be obtained from equations (147) and (148), 
remembering that II, = II, = 0 and IIz = Y : 
Also, the above equations were for a generalized 
cylindrical field, while we are now restricted to 
Smrculariy cylindrical functions and therefore, 4,=4,=], 
Me=f 4» and yer , U*=6 , 272 - Note also 
that, due to the exponential e@ */h2-<ut , aifferen- 
tiation with resvect to Z multiplies the function 
by +i, . For the transverse magnetic field, the 


results are? 


gp x tR 


le 


dh O 
(150) Eg? = +th 


® 


En! = (#-Rp 
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Coe eee 
( a 
(61) He = ab $F 
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fe 7 O 


Similarly for the transverse electric field 





ae uw aw 

rr “rf oo 
(2) | OV 

(162) Fo =-Auw SF 
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(163) Bees 


\ 


He = (AK) 
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When initial conditions are prescribed over given 
olane or cylindrical surfaces, a solution can be 
constructed by superposing these elementary wave 
functions. For fixed values of Ww (or #) and h, 

the field in cylindrical coordinates may be 
represented by the eauations (164) and (165) below, 

in which the q, and 4, coefficients must be determined 
from the initial conditions. 


am 


(164) 


E,> &-#) > GoW 
Mi=a—- F 
2 = 2 
H= 44 Sn Uy + By Un 
Na =-0o 


(165) He-# > anny ~ 2 abn Yo 


Hz= (R° —#) Sie Un. 
IN =- 0 
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Propagation alone a circular cylinder 


Consider 2 circular cylinder of infinite length, 
of radius a, of material with propagation constant 
ki imbedded in an infinite homogeneous medium of 


propagation constant k The field must be finite 


2° 
at the center and conseauently the wave functions 

within the cylinder will be constructed of Bessel 
functions of the first kind. Outside the eee 

the field will extend to infinity and consequently 

the functions must be constructed of Hankel functions. 

In the equations below x *4is substituted for (£°- £2) 

and x 1s substituted for (,-#°) . The function 
eine take raw is abbreviated to Fn. Summations 

are in each case from n =-°%o to n 2+“ nay eis 

prime above a cylinder function indicates differentiation 
with respect to the argument Oge . lhe super- 

scripts 1 or e over the coefficients 4g, and 4, 


indicate points interior and exterior to the cylinder 


respectively. 
Et= > [Saneanas— Sard, b0) o| Fn 
i pt, AUW q & 
(166) a S [aR Sar) an + "% J, Kir) bs \Fn 
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He= > [2M donai+ 451609 bi] 
(167) He ae > | 4 Un (4) ai — 2 a Or) bi| Fr 


= y (In(Ar) by] Fr 


e- La 

bE, = S (sf H,, (a1) as~ aa? He He at) b8 | Fn 
5 a 

(168) ES=— > [26 Hl carat +242 at! Ou) af Fr 


= > {in "(GN ae i 


=> ae, HY (ar) a8 4 hi (ar) ile 


(169) HE = S [8 Ht enat—-SE We ed o6| Fn 
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The coefficients of the expansions end the 
propasation factor h must next be determined. At 
the boundary where ris equal to the constant radius 
of the cylinder, a, the tangential comvonents of the 
field are continuous, which condition supplies a 
relation between the coefficients. From the 
continuity of the tangential components of E, the 
following relations are found, where the substitutions 


have been-made xX,@ =P a 


(170) 


; f Sd | )' 
TE, dy (Pag + 38 (Poe = PH w)ae + eH OBE 


elf) dn(Pae = it (wag 


From the tangential components of H, similar 


relations are found as follows: 


(172) 


8 I te)as ~ %& Tyee = 2% Bdge — 74 Hl eles 


4, WP Ay War 








63. 


IN 


(173) Ty lbp = Hi (was 

The above relations constitute a homoseneous 
system of linear eaueciicns satisfied by the four 
coefficients O;! : ba a Da . There will be 
a& nontrivial solution only if the rank of the matrix 
is less than four, which will be true if its 
determinant is zero, The pronazation factor is 
determined by the condition that the determinant 
of (170) throuzh (173) shall vanish. 

Expansion of this determinant leads to the 


transcendental equation 


(174) 


yf AL ee) fr" a. 
eo - a pte | A al P i Mar i ME (on 7 a) 


et 


The roots of (174) are the allowed values of the 

propagation factor, h, and so determine the 

characteristic or natural modes of propazation. 
Thus, there has been develoned an equation 


whose solution enables caleuletion of natural modes 
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and the propazation constants for long rods. Equation 

(149) expresses the true result of Maxwell in an 

infinitely long rod imbedded .in a perfectly 

homogeneous medium. Equations (157) and (158) are 

valid solutions to this equation. By use of the 

continuity of the tangential components of » and H 

at the boundary of the rod, it was possible to use 

(157) and (158) to find relations which enable one 

to solve simultaneously for the propagation constants. 
Solutions to (174) have been obtained and it 


has been shown?? 2° 


that they verify exoerimental 
results within the accuracy of experimental procedure. 
In the next chanter, a nomogram will be presented 
which provides solutions for the case of dielectric 
rods in air or vacuum, 

It should be emphasized that the preceding 
development assumes that the rod is long enough 
to neglect boundary conditions at the ends, and that 
mae radius, a, is not microscopic as far as the 
external wave is concerned. The question as to what 
size rod is reauired to justify the assumption that 
the radius is not microscopic is most pertinent. 
If the radius is small compared to a wavelength, 


it is probably in error to assume that Hankel functions 
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= 


are required for all waves exterior to the rod 
while Bessel functions of the first kind suffice 
for internal waves. 

Having obtained the provagation constants, 
available ovublished literature indicates that all 
authors proceed to determine radiation characteristics 
on the assumption that the propagation is truly 
constant throughout the rod. The experimental 
evidence Crete this is not true in short rods of 
diameters of a wavelength or less has avparently 
been neglected, and available published theories 
on vowers show serious defects when confronted with 
experimental evidence. For exammle, experimental 


8,12 show that power On @xis 1s linear with 


curves 
exvosed length of dielectric for short antennas, 
The published theories do not permit such a 
ied ton. 

A satisfactory solution to this problem from 
a theoretical standvoint cannot be completed in 
the vresent work. The experimental data presented 


may provide a stepping stone toward its ultimate 


achievenent. 
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CHAPTER 2, 


APPLICATION OF =LECTROMAGNETIC THEORY TO 
CIRCULAR DIZLECTRIC RODS 


In Ghapter 1 an equation was deveiloped 
expressing in circular @oordinates the naturel 
modes of electromagnetic waves in an infinitely long 
cylinder imbedded in an infinite homogeneous medium, 
The next step is to svecialize this to the case of 
@ dielectric cylinder imbedded in air or vacuum. 
Furthermore, particular solutions of (174) are 
desired, and to achieve this, the equations will be 
put in a form more suitable for use with available 
tables of cylinder functions. 

The eauation (174) may be expressed in Modified 
Bessel functions a functions) by use of the 


relation 


(175) ki) =i be ™® Ho (a) 
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(176) H (v)=-iFe aa (-inr) 


and 


ay a 
(77), (w) =-# O° KK, Ge) 


Letting -iv = q, then 


Hp cee 
(178) As Hy (ar) G Knl(Z) F 


We let the above expression be -g, by definition, 
to simplify notation. Also let 


ae is a _ 


Substituting these relations in (174) 


(180) 


(u, fF tu gX # f t & 9) = am AG (<3 + +t) | 


” 
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since we are only concerned with dielectric rod 


and air or free space, we may consider 
(181) 4 =H, =H 


The constants k and h were introduced in our 
elementary harmonic solutions of the scalar wave 
equation in a cylinder, equations (143-145). From 
this, k was equal Vw'4f -iweu , and we may now 
set 6=0 » since the conductivity of dielectric 
and air are negligible. Thus A= wlye, ,and 
KR, = w¥ué, . In the notation introduced with 
equations (166-169), X* =A°-A* and Y=h-K , 
Since x,= = and X= ¥ =4L, it follows that 


(182) 


Rr + £ = wu, 
(183) 
Ae ~ $5 = wu Es 


Where the subscript on é, has now been changed to 


zero to indicate the snecial case of free space. 
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Eliminating h from the above, the following 
is obtained 


ass ESE = wu (6-4) 


at 


Let us now set for the relative dielectric constant 


(195) € = & 


and use the relation 


(186) C= TR 


to put eauation (184) in the form 


Wa . aa _ \;/Pttg> 
oes = oe eae 


which will be one of our principal relationships 
to use in solving (174) for the propagation constants. 
To obtain the other vrincival ecuation, we 
substitute into (180) the value of h in terms of 
the other constants as civen in ecuation (182), 


substitute for k, and Ss their cormesponding values 


- 
as above, and substitute for Ww from (187). After 


simplification, using (186), the result is 
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(188) (Ef +g)(f +g) = a (G4 lft) 


For a givenn,w, and é, the last two ecuations 
can be solved for p and a which define the mode 
Brevazated along the rod, and sonable calculation 
of the other constants. 

From equations (182) and (183) some interesting 
deductions can be drawn concerning the lonzitudinal 
propagation factor, h. ‘We note from (183) that if 
q is zero, that h is equal to wV¥uvé , and therefore 
the wave is traveling with the velocity and wavelength 
characteristic of free space. If p is zero, (182 ) 
shows that h is ecual wifué , and the wave is 
traveling with the wavelength characteristic of the 
dielectric. For a given dielectric, freauency, 
and diameter, equation (1287) shows that the sum 
(p- + gd 2) is a constant. It will later be shown 
rat pis subject to less variation then q. If 
@ is small, which is the case for small diameters, 
the K (aq) becomes very larze compared to J fp). 

Since the wave functions outside the cylinder were 
constructed from a TunetLons, wiien were Later 
converted to K, functions, it is apparent that the 


energy outside the dielectric is large compared to 


_ that inside. Also for small gq, the wavelength and 





ale 


velocity characteristics are near that of free 
svace. Conversely, for large diameters, q is 
large relative to p, and K,,(a) becomes very small 
compared to Jp). Thus the energy is contained 
mostly inside the dielectric, with a relatively 
maya itlm Of enerey on the outside surfaces. The 
wavelength in this case anproaches the characteristic 
of the dielectric. Qualitatively, we mizht expect 
in view of this, unless the dielectric is ideal 
and lossless, that losses increase rapidly with 
increasins diameter of the rod. 

An important factor needed for the design 
of dielectric antennas is the ratio of free svace 
wavelength to the wavelength \propacation, or 
aoparent index of refraction. We have' noted, 


where conductivity is zero, that: 


(a) the vropazetion constant, characteristic 


of the dielectric, k, , is equal wYuUé, 


J 


(b) the propagation constant characteristic 
of free space, 


Ky = Kyo = wlu& = Gy 


(c) the Longitudinal propagation factor for 


the combination is 


peer 
h x 








ie. 


me 1Oollows that 





2 " 
(189) Aes ree A 


where no is used as a symbol for the apparent 
index of refraction. In order to obtain an explicit 
expression for no in terms of the other constants 
in these equations, divide caquation (190) by Ke 
and note that 


eee 
A. = e 


zi 


Substituting for the numerator its quivalent given 
by (188), and extracting the squere root, we have 


n, ieeeerns Of € , p, and q. 


Ce +9)" 


(191) Va = (pe + Vy 





Con 


imeoen is in @ convenient form for obtaining n, 
when p and q have been determined. 

Curves of calculated avyparent indices of 
refraction for various pene tere and dielectric 
constants for the dominant mode (n = 1) are given 
in Figure 5. These are in reasonably close agreement 
with expverimental values obtained by MeKinney?®, 

It can be readily understood that this 
aonvarent index of refraction is of fundamental 
mmporvance in the design of a dielectric antenna. 

It has been shown eae that the optimum length 
moea cylindrical dielectric rod, to maximize 
radiation on axis, is that length which causes the 
wave to emerge from the end delayed by a half period 
behind the externel wave. This length has teen 
shown®s 11,13 to6 be 


(192) 4, = ao 


which is obtained as follows. Consider that the 

time for the wave phase to vass throush the dielectric, 
ta must be greater than the time for the phase to 
sweep through the same distance in air, tg, by one 
half period, - ° 
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Substituting c= “4, rand v, = j6), = w de 


and solving for #,, ecuation (192) is obtained. 

The above expression for the optimum length has 

‘been verified expverimentally by several investi- 
gators!»ttl, 

In the discussion so far, the mode of the 
electromasnetic waves, which is determined by the 
n in the precedins equations, has not been 
specified and the equations hold for any mode. 
it was shown in Chapter 1 that, if the wave is to 
be neriodic in ©, that n must be an integer. For 
2 practical dielectric antenna, it is necessary 
SO specify this further. The mode normally used 
mor Dropesacvion in the wave guide and for exciting 
the dielectric is the so called dominant mode, 
which is tne mode of provagation having the lowest 
cut off frequency. Shelkunoff has aoune by 
allowing q to approach zero and evaluating the 
remaining terms in the equation, that the n = l 
mode has a cut off frequency which is theoretically 


zero. For n zero or greater than 1, it is also 


shown that the cut off is not zsro. Therefore, 
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the dominant mode is n = 1 and this is the only 
mode considered for the remainder of this work. 

The method of solution of equations (17) 
and (188) must be either numerical Ormera vaca. 
and both of these methods are cumbersome and 
tedious. In order to facilitate obtaining ravid 
solutions, at a small sacrifice of precision, 
nomogzrach charts A and B have been vrepared. See 
Figures 6 and 7. 


Chart A, Figure 6, is a nomogrem of the equation 


2 + qé = constant 


-_ 


and Chart B, Figure 7, is a nomogram of 
EF +gXt +9) — pr 4+ 4) =0 


To obtain a solution, we seex values of pv and a 
to satisfy both equations. It is assumed that the 
relative dielectric constant, €, the radius, a, and 
the free spece wavelength, Ao , are known. To use 
moe accompanying charts, proceed as follows: 

Bovcurate tne constant, pe + qe , by Slide rule, 
from the eavation (187). Thus 

ap © egy 
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and locate this noint on the center column of 
Chart A, A straight edze pivoted at this point 
will give the values of p and q allowed for the 
eaven dielectric and ciameder: Assume a value for 
q and vlace the straight edse at this point on the 
left column, keeping it also at the pivot point 

of the center column. The corresponding value of 
yp is read where the straicht edge intersects the 
right column. 

Now look at Chart B, and marx off or visualize 
an € line, corresponding to the given value, running 
straight and parallel to the dash dot lines, Find 
the assumed value for q a@long this line, as indicated 
by the dashed q lines. At this point, thus located, 
read the correspondinz value of p by eye inter- 
~Olation between the solid p lines. Compare tne 
values of p as given by Chart B with that obtained 
from Chart A. If the value of p as given by Chart B 
exceeds that given by Chart A, then the assumed 
value of q was too large. Assume a lower value of 
q and again find the corresponding p values from 
Charts A and B as before. ‘when the two charts 
give identical v's for the same q, the solution 


4s correct. 





: 
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The above is a trial and error method, but 
it is easily verformed and converges very quickly 
wemoae solution which 13 sufficiently accurate for 
many purvoses. It msy hannen that a solution may 
pe required that does not fall on the chart. It 
was imoossible to meke a chart which was readable 
and of reasonable size without considerable 
me-ouriction on the range. However, the range 
selected for the accomzanying cnarts is believed 
to be useable for most of the practical problems 
anticinated at this time. If the required solution 
does not fall on the charts, numerical methods should 
be used. In order to facilitate these, the values 
Sea end ~¢ 42s functions of p and gq, respectively, 
may be obtained from the plots of these functions, 
given in Figures 8 and 9, These plots may also be 
useful to refine the gravhic solutions obtained 
from the Charts A and B, if additional accuracy is 
desired. 

Having solved the enuations for the waves in 
mecrrcular d.electric rod in air, we can now 
readily calculate the optimum length for a given 
G@ielectric constant, diameter , and freouency 


(equation 192). The next vroblem is to determine 
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what gain can be exnected from the antenna. There 
is no e@ccepted method for ce@iculating thercain, 
although several theories have been proposed’ +? 1? ° 
It is for this reason that we turn to experiment 
to obtain the gain to be expected for the given 
eonditions. The next chapter will describe the 
equipment and procedure used to obtain the needed 
experimental data. After describing these 
measurements, the data will be discussed and an 
attempt made to correlate the results with Mersten's 
theorytt, | 

Lt will be noted that the preceding theory 
has dealt exclusively with the circular cross section 
rod antennas. The case of the rectangular cross 
Becton 2S much more difficult due to the fact that 
there are discontinuities on the edges, which are 
cumbersome if not impossible to account for. Further- 
more, there is no convenient coordinate system in 
which to set down the mathematical expressions. 
No published theory of electromagnetic waves in 
‘dielectric rods of rectangular cross section is 
available, and consequently for the design of dielectric 
antennes of these shapes, the exnerimental data 
must necessarily be used. It will later be evident 


that the nerformance of rectangular rods is very 


de 
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definitely different from that of circular rods. 








On 


CHAPTER 3, 
DESCRIPTION OF EQUIPMENT, PROCEDURE, AND @XPZRIMTNTAL DATA 


The objective of the exverimental work was to 
determine the gain on axis for dielectric antennas 
excited in the dominant mode, using optimum length 
rods in each case. These were measured for cylindrical 
rods of varying cross section diameters, and for 
rectangular rods with varying sizes in the E dimension, 
holding the H dimension constant at the width of 
the standard 3 centimeter wave guide (0.9 inches). 

A wavelength (A.) of 3.2 centimeters, eauivalent 

to a frequency of 9.37 x 107 cycles per second, was 
used for the reason that the components were of 
convenient size and the equipment for use at this 
frequency was available, The equipment tsaoshown in 
Puree 10-14. Two dielectric materials were used 
for each antenna: Polystyrene (€=2.56 ) and Teflon 
cae ). 

A standard Kilystron oscillator was used for a 
transmitter, but square wave modulation was imvosed 
on the carrier freouency to facilitate detection in 


the receiver. The transmitter operated from a 
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regulated power supply which was fed from the 60 
cycle 115 volt power lines through a voltage 
regulating transformer. The receiver consisted of 
a crystal detector, amplifier, and voltmeter. As 
a cneck on receiver operation, all measurements 
were repeated usinz a Pickard and Burns Bolometer 
Amplifier, model 100. Results were in substantial 
agreement with those revorted and will not be 
repeated. 

The various sizes of dielectrics were fitted 
into sections of wave guide, made esvecially for 
the purpose, tapered down to standard size and 
fitted with a standard choke. See photographs. The 
part of the dielectric which fits into the suide was 
tapered in each case, in order to minimize reflections 
and reduce standing waves. The elimination or : 
correction for the standing waves proved to be the 
most difficult problem to overcome, The taper 
reduced the voltage standing wave ratio (VSWR) in 
most cases to below 2.0. In every case this was 
carefully measured and a correction applied to obtain 
the actual transmitted power, using tne well known formula 

(193) 


ait ee 
P. “SGA 





yhere P, is the transmitted »over 
Ps is the incident power 
r is the voltage standing wave ratio (VSWR) 


several methods were used to make the measurements. 
To begin with, the dielectric was used as the 
GEansemiovins antenna, sinmee recSiver power was 
inadecuate to accurately measure the VSWR, The VSR 
was measured at the transmitter using a standard 
standing wave indicator. Later it was found more 
precise to use a calibrated attenuator together 
with a slotted section, crystal detector, amolifier, 
and voltmeter, attenuating as necessary to equate 
the voltage from the hich voltage point to the low 
voltage point. The method by which the final 
measurements were made used this equipment, emnloying 
an attenuator which had been accurately calibrated 
by the National Bureau of Standards. 

As a further check on the accuracy of the 
measured gains, the various dieleetric antennas 
were each vut et the receiver with a svecial matching 
section, and the measurements reneated. For each 
reading, the matching section was adjusted to 


maximize readings by the reduction of the standing waves. 
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Results for the dielectric antennas un to one 
wavelength in diameter or in the # dimension were 
in agreement with those measured by the vorevious 
two methods, Gains measured for larger antennas 
mere somewhat less, which e@m be accounted for by 
the fact that standing waves still existed to an 
appreciable amount cespite the matching section, 

Some slisht difficulty was encountered with 
fer tcetions within the room. these were reduced 
as much as vossible by the strategic use of an 
absorbing screen. As a check on errors which might 
be introduced due to reflections, the distance from 
the transmitter to receiver, normally of the order 
of one hundred wavelengths, was changed in quarter 
wavelength Increments, noting nower detected at 
eacn range. If anoreciatle reflections were 
interfering, there snould be alternate addition 
and cancellation of power observed. It was found 
that detected nower was not anovreciaoly different 
for these vositions 4 quarter wavelength anart, 
and it wes concluded that any errors from reflections 
were of little significance. 

Before meesurins the gains it was necessary 


first temenca the optimum lengtn of the dielectric. 





This was done experimentelly by observing relative 
power transmitted, first with the dielectric almost 
flush with the mouth of the guide, and then out 2 
inch and so on with longer lengths using increments 
of about = Ro . Near the veak of the transmission 
smaller increments were used to obtain the exact 
optimum length. Having the relative power trans- 
mitted at optimum length, the test antenna was then 
immediately replaced by a precalibrated standard 
horn and the vowers compared. This enavdled 
calculation of the relative gain of the test antenna 
with respect to the standard horn. 

The gain of the standard was carefully 
measured in order to be able to determine the 
absolute gain; that is, the ratio of the neak intensity 
of the given antenna to the intensity of an isotronvic 
radiator when both are radiating the same total 
pover. In terms of the effective area, Ae, the 


gain, G, has been shown” °°? to ve 


(194) a 


@ 


The effective area is the ratio of the vover 
received at the antenna to the nower ner unit erea 


in the incident wave, 


SB. 


The absolute gain of the standard horn was 


measured by several of the methods which have been 
Dipeeeo 20 


s 


develoved and ynuolished - One of these will 
be briefly described here, which makes use of tvo 
identical antennas, The standard horn was attached 

to the transmitter and an identical one used on the 
receiver. Under these conditions it is easily showm@ 


that 





2 
(195) i = (= 

where fF, is the received power, fF, is the transmitted 
Sewer, and d is the distance between the transmitting 
and receiving antennas. To determine the ratio of 
receivec »vower to total transmitted rower, the 
transmitter and receiver were couplea directly 
together, except for a calibrated attenuator between. 
With a large but accurately measured amount of 
attenuation inserted, a varticular receiver voltage 
was noted, Then the two antennas were attached 
respectively to the transmitter and receiver and 
senarated by an accurately measured distance 

(several feet), and the attenuation of the transmitter 


reduced until the same receiver voltage was obtained. 
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The difference in the amount of attenuation in the 


two cases gave the ratio of PL to Pr 1gl (6iexe)a. \efedLS) 


R 
from which G could easily be calculated, 


Presentation of Data 


Three alternatives in the vresentation of 
experimental data have been considered, and these 
will be briefly mentioned to show the reasons for 
the selection of the one given. 

(a) Renort the relative sain of the dielectric 
antenna compared to the gain of the feeding wave 
eulcge with the dielectric removed. This would 
permit a direct comparison of the nower on axis 
with and without the dielectric which would show 
clearly the effect of the diclectric in directing 
the beam. However this method was excluded, because 
the experimental data would not be reproducible 
unless the feed structures were identical. Gain 
of oven mouth feeds are not vrecisely predictable 
and there would be no way to obtain absolute 
gain or sfifective area from these data, 

(b) Report the absolute gain of the dielectric 


antenna; that is the cower on axis relative to an 


Cn: 


isotropic oscillator emitting the same total power, 
This is better, but here asain variations in the 
metallic feed structures would cause variation in 
the absolute gain. Furthermore, as shown by 
equation (194) above, the scales would have an 
ertitricial multiplier of 47. 

(c) Revort the effective cross section area 
relative to the ohysical cross section area. This 
is the method chosen because: (1) this presentation 
eliminates the undesired effect of varzation in 
feed structures; (2) as cross section areas are 
ablowed to grow, the effective cross section areas 
also srow, and these tend to become equal. In other 
words, the captured energy tends toward that 
physically intercested by tnis area. Therefore, 

a plot of effective area versus physical area will 

verge toward a straight line with unit slope passing 

throush the origin regardless of the dielectric 

used; and (3) an extremely useful piece of information 

concerning any antenna is its effective cross 

Peaemwon area, From this figure it is a simple and 

direct calculation to obtain tne total energy 

gathered by the antenna in a field of a siven strength. 
This third method of presentation, in addition 


LO DYOViGiwew= tre engincer the very useful information 
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of maximum capture area for various dielectrics 
of various sizes, illustrates the variation of 
capture area from the actual area in the sizes 

of the order of a wavelength, and shows how these 
should converge to the 45° line as @imensions 


increase toward ontical sizes. 


Choice of Scales 


Consideration was given to reporting dimensions 
Girectly in centimeters or inches. This would 
enable an engineer, employing wavelengths the same 
as used here (3.2 cm) to use the data directly. 
However, this was rejected because the value of 
these cata is not limited to any particular wavelength. 
Dimensions are therefore sziven in non dimensional 
form; that is, lengths in wavelengths and cross 
section areas in square wavelengths. Although this 
reouires conversion of practical units to this system 
to enable usage of the data, it has the advantage 
that the data are apolicable to a very wide band 
of wavelengths. 


A plot of effective area versus actual cross 
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section area for the cylindrical dislectric antennas 
measured 1s siven in Figure 15. Figure 16 is the 
same plot for rectanzular cross section antennas. 
Figures 17 and 18 are also enclosed for completeness 
giving the absolute gains of cylinders and 
rectansular cross sections, resvectively, versus 
egos Sectlon area. 

It will be noted that measurements were not 
Sovained for cylinders of larze enougn diameter 
to shor convergence of the effective area toward 
the actual area. In the larger diameter circular 
cylinders, the mode tended to shift to the m=0 
moae whoienm gives a null on axis. No way to prevent 
this was devised. There was also a practical limit 
to the measurements on the low side, with the time 
and facilities available. For small cross section 
rods thé optimum length becomes unmanaseably longs, 
and the thin rods tend to sas. Furthermore, these 
thin longs rods reauire a long range in order to 
nut the field meter in the far zone. 

The present techniaue does not permit extension 
of the range of measurement for circular rods, 
and the vossibility of existence of TM modes in 


these dielectrics is a disturbing factor. The 
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aovarent irregularity of the polystyrene rods may 
be caused by this nhenomenon. In order to reduce 
the probability of spurious modes in the exciting 
madiaver for the larger areas, conical horns were 
also tested with polystyrene cylinders, with the 
following results 





Horm Antenna 

Horn diameter Physical Area Canture Area 
ei" 2.5 22 3.78 A: 
3i" 5.25A; 5.24 AS 


Here the trend of antenna effective area toward 
actual area is noted. The horns were far from 
optimum from the voint of view of radiation pattern, 
however, and the lens antennas were also probably 
not ovtimum. Therefore, convergence of opvtimum 
cvlinders toward the exvected asymntote would not 
be closely realized until the experiment reached 
much larger areas than the above figures indicate, 
In the rectensular rods, no tendency of the 
mode to shift was observed, but other nractical 
difficulties prevented extension of the rance of 
sizes. In order to prevent cut off in the antenna 


feed, when using small sizes, the H dimension was 
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set at a constant value of .71,A.. Small areas, 

then, would have reauired very lons and thin sections 
which would have caused the same difficulties mentioned 
above for circular rods, The larger cross section 
areas, where the E dimension was large compared 

to the constant H @imension, gave radiation natterns 
which were relatively narrow in the plane parallel 
to the © dimension, and relatively broad in the 
vlane »varallel to the H dimension, thus giving 
reduced gain on axis. It was found that if the 
cross section area had been made more nearly sauere, 
a radiation pattern more nearly circular would heve 
been produced, and thereby increasing gain on 

axis. The gain on axis, being a function of the 
radiation pattern, then, is denendent on the shape 
of the cross section. For small cross sections this 
effect 1s small, but for cross sections of a weave-~ 
length or larger the effect is very pronounced. 

A way to avoid this would be to proverly shane 

the front surface of the dielectric, or to taver 

the rod, Neither of these were attemvnted here. 

The end surfaces of all the antennas tested were 
vlane, 


it was evident that the larser rectangular 
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dielectrics would give hisher povers on axis if the 
cross section was anproximately square. Therefore, 
as the wavelength region is approached in Figure 16, 
the larger vhysical cross section antennas apnear 
at a disadvantage. Perhaps a better comparison 

in this region would have been to compare the antenna 
capture area to the capture area measured for the 
same metallic feed with the dielectric flush with 
the metellic anerture (zero length antenna). This 
method of revorting would have brought the larger 
areas more in accord vith the following table, 


measured only for pnolystyrene. 








Feed Horn Physical Area Canture Canture Area 
. Ares 
Dimensions a Ac of of Antenna 
Feed Horn 








ROG" 1.02 0.64 1.03 
i; x 13" Tee 140 ey 
eee 3c 4.76 4.75 5.3 


This table shows the trend, predicted 
oreviously, of the capture area of the ontimum 


length antenna to avproach the vhysical cross 











section as the latter becomes many wavelengths in 
each direction. The shaping of these larger 


antennas was only attemnted avvroximately and 


somewhat larger capture areas would be expected 


if ontimum shaves had been employed. 
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CHAPTER 4, 


DISCUSSION OF RESULTS AND ATTEMPTED CORRELATION WITH 
EXISTING THSORIES 


The curves of ieee area versus cross section 
area, Figures 15 and 16 of Chapter 3, show that 
Comsidenable increases in effective areas over 
actual areas are obtainable if the physical areas 
are small (less than a square wavelength). These 
effective areas are auite high for small physical 
areas and reduce toward the vhysical areas as sizes 
become larger, The trend is avparent for the 
rectangular dielectrics, but could not be demonstrated 
with circular cross section because of shifting 
of the whode within the dielectric despite excitation 
by the h = 1 mode, 

As previously pointed out, the design of the 
dielectric antenna above a wavelength in cross 
section area is still not understood, The gain on 
axis, on which the exnerimental results are based, 


Ls a function of the snape of the total nattern. 
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The gain is measured by samplins this pattern on 
the prolonged axis of the dielectric. This svread 
in vattern is particularly avparent in these data 
in the sizes immediately above a square wavelength. 
In this region ray optics play an imvoortant role; 
and yet customary treatments by ray optics fail 
because of the considerable vohase gradients due 

to the adjacent metallic structures supporting the 
radiator. In fact, such a treatment would involve 
serious difficulties. In the vresent project, the 
prooviem would be further comnlicated by the absence 
of an accurate means to determine the comolste wave 
characteristics of each or tne various feed horns 
used. Proper lens shaning would necessarily have 
to be based on a knowledze of the chase front of 
the horn to be lens:d. However, to show the trend 
beyond the one wavelength region toward the exnected 
45° asymptote a few measurements were attemvoted 

in which the face of the dielectric was shaned by 
means or trying many combinations of various si.e 
blocks so as to obtain a stenved vorofile., The 

best gains were found in these larzer shaves when 
the profiles resembled ovtical lenses in physical 


shane, and the larger the area the more nearly the 
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true optical shane is reauired to optimize gain. 
Tois is in contrast to smaller cross section areas, 
where shaping of the frome tace Had little eitect, 
and exverimental measurements were very consistent. 
Evidently the vnhase gradients within the smaller 
dielectrics play a much less significant role. 
Because of the fact that these data are based 
on unoptimized vrofiles, it is apvarent that gains of 
Optimized shanes would be higher than haere reported, 
and convergence to the 45° asymptote would occur much 
further out. It is exvected that when better theories 
miro available for the design of lenses in the wavelengtin 
region, curves similar to those of Figure 16 could be 
exnected to converge monotonically to a line of unit 
slone, and to reach the asymptote only after dimensions 
are such that the areas are many sauare wavelengths. 
One exveriment was attemnted some years aso to 
"lens" a 30 inch parabolic averture on X band, and it 
was not found nossible to detect any increase in 
effective area, As the anerture grows, ontical lenses 
are obtained, where effective area is strictly provortional 


5 


to actual area, Ina recent report of exneriments 
sponsored by the U.S. Army Signal Corps, tests were 
revealed on a spherical Lucite lens of cross 


section area 985 square wavelengths. Gains 
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reported, after correction for the losses occurring 
in Lucite, indicate an effective area of 635 square 
wavelengths. The aveawene discrepancy is due to 

the fairly narrow illumination used in this 
experiment. The illumination duplicated the usual 
radar dish illuminations and resulted in the 

expected averture utilization factor of 65% 

commonly employed in radar dish antennas. This shows 
that effective areas actually do converge to the 
moysical areas utilized in large sizes. 

Converging lenses have been made for the 
wavelengths considered herein. While these appear 
to give a tremendous power at their designed focal 
point, a more thorough discussion would show that 
they merely operate as converging devices, obeying 


to a first approximation the half power law 


4 = power angle of radiation 
0 


They do not increase the effective area of the 


aperture. 


VMarsten's Theory 


A very excellent theoretical treatment of the 
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aot erate antenna of circular cross section is 
given by Marsten in a recent revort!+ nubli shed 

by the U.S. Naval Research Laboratory. A circular 
dielectric rod is excited in any sinsle mode, and 
the ecuations used for the waves 2re equivalent to 
those given in Chanters 1 and 2. rom these are 
calculated eouivalent electric and magnetic current 
distributions over the surfaces, in accordance with 
Schelkunoff's Eauivalence Theorense”. Then the 
radiation field is computed using the method of 


echelkunoff, first for the field arising from the 


eavuivalent sources at the sices of the rod, ea 





and then for the equivalent sources at the end of 
the rod, £% 7 Certain factors common to both 
expressions are dropoed, since only the relative 
field is » Of interest, and the final expressions 


for the © components are given as follows 








EY = ~cosng 2 2% 79 (y scose+ld,,(1) — E82 J (L) + 


mM 


(196) 
(,scose-I)J,,,(L) | 








oo. 


(196) continued | 
ee ee je * Ph . +u4 8) J, (L)-J,(e) + 


er ai P iL I(t) 4, ~ pd, (ICL) 


The expressions for the other comronents of the 





field will not be reveated, since only the power 
on axis is of interest here. 
The notation is the same as given herein, 


exceovt as follows: 


J ee = Vay 
: 
f= we = 4) 
fp ------- = length of rod. Since optimum length (2) 


; Sik 
4s here considered, substitute 2m-l) 
ve 
ml, Cos @ - 4) 


a 
it 


On axis Cos@ = 1 and for ontimum length X= - 


= $a sno On axis [|=0 
A and B are functions which are always finite 


and non-zero 








LOS 





— BH In (IK @ ~ B* Kingda) 


a a re sem ee ee ee 


mh (p+ ze) Sn (P) Kn (G) 


ay 
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which upon simplification and substitution from 


equations (187) and (194) becomes in our notation 


lh 
~ 1 &@f+p _ ,§ V 
oT (+9) K 


By actually substituting numbers, it is found 
that §& has a nlus value for all dimensions under 
consideration, and so the vositive sizsn snould be 


used after extractins square root, 


Assuming n=l, 6-0, y = A , ana p=0 
the above equations reduce to 
(197) 
o Se eee 
6 = Tr (ined i 
(198) 


) 
cS == © (gay a UAL Vimgt!) 16 tg LAs 
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A plot of effective area for optimum length comouted 
as above is given in Figure 20, where the abscissa 
is taken as tne physical area, 

The gains calculated by Marsten's theory show 
the same general trend as the experimental data 
of Chapter 3 show. The vrincinal excentions are, 
first, that all are much h¥aher in value than could 
be found exnerimentally, and, second, it is noted 
that the theory does not hold at all when dimensions 
are small enough so that n, approaches 1. HZaquation 
(202), the expression for the contribution from the 
sides, has a demominator, [(t.-!) which approaches 
zero as n approaches 1 with decreasing diameters, 
The term 7,(rF) in the numerator also approaches zero 
as the diameter goes to zero, but this occurs much 
more slowly. The result is that vowers on axis as 
indicated become extremely high for sm@ll rods, These 
enormous powers do not conform to experimental 
observations. 

The question of identical phases for the 
radiation from the side walls and from the ends 
does not appear to be correct from a physical viewpoint. 
The phase within an optimum length rod emerges one 


half period later than the wave travelling the same 
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distance in free snace. Along the rod the vhase 
relative to the end must certainly lead by something 
less than a half period but more than zero. If, 
in the theory, some adjustment to these nhases 
were made, it is possible that the theory might 
fit the experimental findings more closely. 

Another interesting point is noted from 
eouations (196) and (197). As the diameter becomes 


very large p approaches zero. ( P=aV/#-% 


and at large diameters h approaches k,). Ifn <=l, 
then J, (p) = 0, and the right sides of both equations 
vanish. We must conclude from this, if Atheory 

holds for large sizes, that no radiation could occur 
from an antenna of very large diameter if the 


optimum length 1s determined as calculated above. 


Wilkes’ Theory 


An exolanation for the peculiar behavior of 
gains in the region of areas of a fraction of a 
square wavelength is siven in a theory proposed by 
vaageeee on the Interaction of Radiant Energy. It 


has been shown experimentally that the effective 
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areas in this region are considerably greater than 
the physical cross section areas, althoush powers 

are definitely limited. Presently published theories 
on dielectric antennas of circular cross section 
would allow an unlimited growth of capture area 

with diminishing cross section areas. 

Wilkes, however, has shown experimentally that 
the power srowth is in reality not unlimited, but 
that the canture area for small cross sections is 
limited. A simplified qualitative outline of his 
aoproach is given below to indicate the method by 
which he arrives at a solution of this hitherto 
unaccountable vnhenomenon. 

As the propagation down a very thin dielectric 
rod is a plane wave, the phase and group or photon 
velocities are assumed by him to be identical. By 
an ingenious techniaue, he has been able to measure 
the ovhase velocities in extremely short lengths of 
dielectric and to comnute the velocities down to 
a zero length. In these thin dielectrics the 
energy starts at its critical velocity, = for 
very short lengths and approaches c for increasing 
lengths. Knowing this rate of change of velocity 


it is possible to conceive of the vhotons rollowing 
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describable ballistic trajectories, This, however, 
reauires the concent of an applied force to vhotons 
of finite mass which exhibit proverties of inertia, 
Since the vhoton mass is at present unknown, except 
for de Broglie's controversial finciness” direct 
use of ballistic pvrincipvles to calculate this 

foree is not pvossible. Nevertheless, with this 
conceot in mind, a consideration of the nature 

of the forces brought to bear on the vhoton reveal 
a behavior closely compatible with exverimental 
observations. 

Postulating the existence of forces between 
the internal and external fields, which we have 
just noted, it follows that there is a transfer 
of work or energy from one of these fields to the 
other. The medium of energy exchange is photons 
of a freauency determined by the excitation. 

Wilkes has found experimental evidence of 
the existence of these forces; in fact his 
calculations and experiments supvort the existence 
of the force as comnosed of two varts. One part 
is a pressure action proportional to the difference 
in velocities of the two fields. The other is a 


function of the total phase lag and appears 
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gravitational in nature. The magnitude of this 
action is found ecual to the radiant pressure and 
parallels in this respect Mandelker's theory! of 
gravitation. 

From Quantum Electrodynamics, Wilkes finds 
that the quanta in the accelerated dielectric 
fields become negative and vresent holes. He then 
allows these holes to be filled in to the limit 
allowed by the statistical laws of a degenerate 
gas, and by means of these caiculations arrives 
at the two saturation vroints shown on Figure 21 as 
the lowest cross sections revorted. these points 
have heen confirmed by his experiments. The other 
voints of Figure 21, namely those for cross sections 
Preater than O.2 square wavelengths, are @ 
repetition of the data obtained in Chanter 3. 
Obviously, if the physical area drops to zero, the 
capture area will also go to zero, and this 
seeminsly trivial statement vermits the completion 


of the curves of Figure 21. 
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Conclusion 


It is believed that no theory based purely on 
Maxwell's eauations will suffice to explain the 
peculiar behavior of gains in the very small cross 
section region. As stated early in this pvaper, 
Maxwell's esuations hold with certainty only for 
large scale electromagnetic phenomena, and when 
dimensions are reduced to fractions of a wavelength, 
they cannot be relied unon., The true exvlanation of 
very small dielectric antennas must come from quantum 
electrodynamics, as proposed by Wilkes, and the 
data obtained in Chanter 3 appear to be in close 
agreement with this theory. 

As dimensions grow larger, effective areas 
of dielectric antennas, if shaped for optimum gain, 
remain greater than the vhysical cross sections, but 
will tend monotonically toward these areas, After 
areas have reached several square wavelengths, 
effective area of optimum shaped dielectrics can be 


closely anproximated by the actual vohysical area, 
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GLOSSARY OF SYHBOLS 


Slectric Field Intensity vector 
Magnetic Field Intensity vector 
Electric Displacement vector 
Magnetic Induction vector 
Current Density vector 


Charge density at any ooint. Also used as a 
general argument of a Sessel function 


Zlecersce DOlari zation Vecuor 
Haeneulc DOlar’ 2Zati on vector 


Slectric Inductive capacity of free space 
defined by D= &F 


Masnetic Inductive cavacity of free space 
defined by B=4,H 


Magnetic susceptibility M=XmH 
Hlectric susceptibility f = fev 


mRolotewe a,c Lecurle Conevany = 
Permeability i 

Vector Potential of the Field 

Scalar Potential of the Field 

Hertz meerer electric Pola Teed olentsal 
Hertz vector, Magnetic Polarization Potential 


Conductivity 
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puUrtace dencmug or ClareS (Or Cuarte "per UnlG area 
Conduction current — 

pauare root of minus one 
Longitudinal propagation factor = Af 
Time 

Base or natural logarithms 


Bessel Function of order n 


An argument of a Bessel function 


Z_(§) Circular eyiinder function, OF Beesel function, 
m 


NCS) 


of order n 


Bessel function of the second kind 


() 
H,. (?) Hankel function of the first kind 


K ni $) Modified Bessel function 


k 


Intrinsic hase constant = Valu +iuew 
Angular velocity = frecuency/zr 


P2OLOpae blo COnsuanl £Or Glelecemic ian. une 
ease 67 4 circular cylinder tna = Ve2-p 


A vropagation constant for air in the case of 
oe Cuseuulatwecy!  ndaer am ais VR R? 


A propagation constant ror dielectric = X,4@ 


BepcOparauLen constent tor air in the case of 
pec tucivar cylinder anvaar = X%.¢ 


Radius of a circular cylinder of dielectric 


A propagation constant for alr = -4w = -ix a = -daVR 





ho 
A 


Xo 


Ors 


Apparent index of refraction for dielectric rod 


{ 
Abbreviation for dn? 
Pda (P) 


I 
Abbreviation for —HKxlg) 
G kn) 


Free space wavelength 
Wavelength of propagation in a dielectric antenna 
Velocity of propagation in free space 


Optimum length for dielectric to maximize 
on-axis radiation 
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